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We report calculations of cosmic-ray proton, nuclei, antiproton, electron and positron energy 
spectra within a “two-halo model” of diffusive transport. The two halos represent a simple, physically 
consistent generalization of the standard diffusion models, which assume a unique type of diffusion 
for cosmic rays in the whole Galactic halo. We believe instead that cosmic rays may experience a 
smaller energy dependence of diffusion when they are in proximity of the Galactic disk. Our scenario 
is supported by recent observations of cosmic ray protons, nuclei, anisotropy, and gamma-rays. We 
predict remarkably hard antiparticle spectra at high energy. In particular, at i? > 10 GeV, the 
antiproton/proton ratio is expected to flatten, while the positron fraction is found to increase with 
energy. We discuss the implications for cosmic-ray physics and dark matter searches via antimatter. 
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Introduction — The observed spectrum of Galactic 
cosmic rays (CRs) is believed to arise from a combination 
of two basic plasma astrophysics phenomena: diffusive 
shock acceleration occurring in supernova remnants, and 
diffusive transport off magnetic turbulence [1-4]. The 
paradigm of CR transport as a particle diffusion pro¬ 
cess has been successful for decades in explaining the ob¬ 
served high-degree of isotropy in CR arrival directions, 
and the abundance of secondary Li-Be-B elements aris¬ 
ing from CR collisions with the interstellar matter. The 
so-called standard models based on this understanding 
employ several simplifying assumptions such as homo¬ 
geneity, isotropy, stationarity, and linearity. The tradi¬ 
tional description is that CRs are confined in a cylindrical 
halo that encompasses the Galactic plane, after being ac¬ 
celerated by supernova remnants to power-law spectra in 
rigidity TZ (momentum/charge ratio), Q oc TZ~‘'. The 
diffusion coefficient is assumed to have a rigidity depen¬ 
dence K oc TZ^ in the whole halo. The combination of ac¬ 
celeration and transport leads to power-law equilibrium 
spectra (j)p ~ Q/K oc TZ~^~'^ for primary components 
{e.g., protons. He, C, 0, Fe), while secondary to primary 
ratios {e.g., B/C, Li/0, Ti/Fe) scale as 4>s/4‘p TZ~^. 
The data on primary nuclei and on the B/C ratio con¬ 
strain the parameters 5 Ri 0.3-0.7 and v Ri 2-2.4. 

New generation experiments, and in particular the Al¬ 
pha Magnetic Spectrometer (AMS), are now probing the 
fine structure of CR phenomenology. The rapid rise in 
the positron fraction e+/(e“-F e+) at E ~10-300 GeV 
of energy [5, 6], in contrast to the standard expectations 
~ E~^, may suggest the existence of extra sources of 
high-energy e^, including dark matter particle annihila¬ 
tion/decay or known astrophysical sources, but these sce¬ 
narios demand a reliable understanding the secondary 
production and their propagation [7-12]. The remarkable 
spectral hardening observed in CR protons and nuclei at 
TZ ^300GV [13-16] may be a signature of new astro- 
physical phenomena, occurring either in acceleration or 
in propagation, that are unaccounted by standard models 


[17-21]. Revisiting propagation, as distinct from acceler¬ 
ation, has also interesting implications for the so-called 
“anisotropy problem” [22]. The anisotropy observations 
suggest a rather small rigidity dependence for the diffu¬ 
sion at multi-TeV energies {6 RiO.15), which is in con¬ 
trast with standard model extrapolations based on the 
B/C ratio {6 «0.3-0.7). Yet a change of diffusion may 
also be hinted by a high-energy flattening of secondary 
to primary ratios [19, 22, 23]. Very recently, new p/p 
data presented by AMS have generated widespread in¬ 
terest [24-28], and some authors have pointed out that 
the interpretation of these data in terms of secondary 
production demands a slower diffusion of cosmic rays at 
high energy [29-31]. We also stress that understanding 
diffusion is crucial, in dark matter searches, for modeling 
both background and signal of CR antimatter. 

In this paper, we report calculations of CR proton, 
nuclei, antiproton, electron and positron energy spectra 
within a numerical implementation of a two-halo scenario 
of diffusive transport. The idea of having two halos with 
different diffusion properties represents the simplest but 
physically consistent generalization of the standard diffu¬ 
sion models. Standard models assume a unique diffusion 
regime in the whole propagation region, i.e., a factorized 
rigidity dependence for the function K. The increased 
turbulence in the disk, compared with that at higher lat¬ 
itudes, leads us to believe that there is a spatial change 
in the CR diffusion properties. This idea was proposed 
in connection with the detection of a break in the proton 
and He spectra [18], but the model was unable to describe 
the sharp structures reported by the PAMELA experi¬ 
ment. The AMS Collaboration has now measured the 
detailed variations of the proton and He fluxes at GeV- 
TeV energies, showing a significantly smoother spectral 
hardening. In particular, the proton spectral index is 
found to progressively increase a.t E > 100 GeV [16]. As 
we will show, a two-halo scenario of CR diffusion is sup¬ 
ported by these data, and it allows us to understand sev¬ 
eral puzzling properties of the CR spectrum. Through- 
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out this paper, we discuss important consequences for 
the physical observables that are being investigated by 
AMS, including implications for dark matter searches via 
antimatter. 

Numerical Implementation — The CR diffusion is de¬ 
scribed by a 2D transport equation with boundary con¬ 
ditions in a cylindrical region for all CR species: 

dt^ = Q + y + ( 1 ) 

where ^ = ip{E,r^z) is the particle number density as 
a function of energy and space coordinates, T = Pena 
is the destruction rate for collisions off gas nuclei, with 
density n, at velocity Pc and cross section a. The source 
term Q is split into a primary term, Qpri, and a sec¬ 
ondary production term Qsec = ; from spalla¬ 

tion of heavier j-type nuclei with rate T^^. The term 
E = —^ describes ionization and Coulomb losses, as 
well as radiative cooling of CR leptons. To numerically 
solve Eq. 1 in steady state conditions {dt'tp=0), we em¬ 
ploy the DRAGON package [32]. We introduced a modifi¬ 
cation of the finite-differencing scheme and of the spatial 
grid of the code in order to allow for a diffusion coeffi¬ 
cient of the form K(TZ,z) = Kq (TZ/TZq)^^^\ Note that 
our function K is non-separable in rigidity and space co¬ 
ordinates. The full propagation region has half-height 
L =4kpc and radius Cmax = 20kpc. We split the cylin¬ 
der into two z-symmetric propagation zones. We call 
inner halo the region which surrounds the disk for a few 
hundred pc (jzj < k =500pc), where the turbulence is 
presumably injected by SN bubbles, while the outer halo 
represents a wider region < \z\ < L) possibly char¬ 
acterized by a CR-driven turbulence spectrum [34] . The 
spatial dependence of S{z) is expressed by <5 = =0.15 

in the inner halo and S = So =0.75 in the outer halo, 
where Kq =10^®cm^s, and TZo =0.25GV set the diffu¬ 
sion in the two zones. The source spectra are of the type 
QP’"' oc 77.“^', with ^=2.32 for nuclei, while proton spectra 
are steeper by 0.08 [33]. For illustrative purposes, we also 
set up a standard model of Iroshnikov-Kraichnan type dif¬ 
fusion, with (5 =0.5 everywhere and source spectra with 
V =2.3, to match the data at E ^ 1-200 GeV/nucleon. 
The local interstellar flux \s p = computed at the 

position 20=0 and rQ =8.3 kpc. In the following, the par¬ 
ticle fluxes at Earth from the two models are presented. 
The solar modulation is described under the force-field 
approximation [35], using modulation parameter values 
4) =550 MV for protons and leptons (relevant for the 
AMS data) and 4> =350 MV for nuclei (relevant for the 
PAMELA data) [36, 37]. 

Protons — The calculated proton spectrum is shown 
in Fig. 1 in comparison with the new AMS data. At 
A > 10 GeV, GR protons are not very sensitive to en¬ 
ergy losses or destruction. Our numerical results nicely 
confirm the trends obtained in our analytical derivations 
[18]. Neglecting interactions, one can see that the naive 
expectation, pp ~ with O = K/L, has to be 

modified using O = {h/ki -|- lo/ko), which depends on the 
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FIG. 1. Top: proton energy spectrum multiplied by E^ ’^. 
Bottom: rigidity dependence of the proton spectral index. 
The model calculations are shown in comparison with the 
data [13-16]. The modulation potential is 4? = 550MV. 


diffusion coefficients fci/o oc and on the half-sizes of 

the two halos l^/o- The GR propagation in both regions, 
with systematically slower diffusion in the inner halo, is 
at the origin of the spectral hardening. The high-energy 
propagation becomes progressively dominated by the in¬ 
ner halo. In the bottom panel, we plot the spectral index 
7 ( 77 .) = d[\og{pp)]/d[log{TZ)]. Our model describes very 
well the smooth evolution of the spectrum measured by 
AMS at 77, ~ 10-1000 GV We also note that, in this rigid¬ 
ity region, the spectrum is never power-law. 

Nuclei — The transport of primary nuclei {e.g., He, C, 
0 or Fe) is similar to that of protons, hence we expect 
a spectral hardening at the same rigidity as for proton. 
A minor difference comes from destructive interactions, 
increasing with increasing masses, that slightly harden 
the spectra below ^ lOOGeV/nucleon. The carbon spec¬ 
trum is shown in Fig. 2. At E GeV/n, its approximate 
behavior is p ^ Q/ {fl hT), where hT is the interaction 
term for a surface gas density hn = 200pc x 1 cm“^, act¬ 
ing in competition with the diffusive term fl. Secondary 
nuclei, like boron in the figure, are expected to have 
steeper spectra and to experience a stronger upturn in 
comparison to their progenitors C-N-0. As a consequence 
the B/C ratio is found to progressively flatten, as shown 
in the figure, to an asymptotic multi-TeV dependence 
^ E~^\ At lower energies, the B/C ratio from our 
model recovers the standard model behavior. A hint 
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FIG. 2. Top to bottom: Energy spectra of B and C, B/C 
ratio, and Li spectrum. The spectra are multiplied by 
The model calculations are shown in comparison with the 
data [15, 38-44]. The modulation potential is 4? = 350MV. 


of flattening in the B/C data is noted in severai stud¬ 
ies [11, 19, 22, 23, 45]. The iithium flux, shown in Fig 2, 
is particuiariy sensitive to CR propagation. Its produc¬ 
tion depends not oniy on C-N-0 coliisions, but aiso on 
tertiary interactions B—s^Li and Be—i>Li. In spite of the 
scarcity of present data, the Li flux is one of the best 
observabies for AMS in terms of detection performance, 
and it is currentiy being measured to a ~ 1% precision 
up to TeV/nucieon energies [24]. These data wiii be very 
vaiuabie for testing our scenario. 

Antiprotons — In contrast to Li-Be-B nuclei, that are 
ejected with the same kinetic energy per nucleon as their 
progenitors, the antiprotons emitted in p-p or p-He colli¬ 
sions have broad energy distributions and large inelastic¬ 
ity factors. This type of kinematics reshapes significantly 
the antiproton production spectrum Their subse¬ 

quent propagation is similar to that of protons, except for 
the presence of the so-called tertiary p-p processes. In 
Fig. 3 the p/p ratio is shown for the two considered mod¬ 
els. In the standard model, it decreases smoothly above 
~10GeV. In our two-halo model, the ratio has a sub¬ 
stantial flattening at E’ ~ 10 GeV - a few 100 GeV. With 
the present data, it appears that no room is left for ex¬ 
otic antiproton components. Recent AMS data have trig¬ 
gered exciting speculations about a possible dark matter 



FIG. 3. The p/p ratio as function of kinetic energy. The 
model calculations are shown in comparison with the data 
[46-51]. The modulation potential is 4? = 350 MV. 


induced antiproton excess [25-27]. A flattening of the 
p/p ratio is an unavoidable feature of our scenario and 
it should be considered when setting upper limits to the 
search for new-physics signals. 

Anisotropy — The GR anisotropy has received renewed 
attention in recent years [22, 52, 53]. From observations, 
the dipole anisotropy in the TeV band remains nearly 
constant to rj ^10“^. Standard model predictions ex¬ 
ceed the observational limits by one order of magnitude. 
To match the data, a diffusion coefficient with i5 Ri 0.15 is 
required. This is indeed the case for the model considered 
here, which gives rj <2x10“^ at 1-100TeV. These esti¬ 
mates depend on the standard formula for anisotropy, 
r] = within the diffusion approximation. We 

also stress that accounting for the discreteness of the CR 
sources may induce large deviations from these estimates. 

Connections with j-ray physics — Although we have 
focused on local observables, our model predicts also a 
CR spectral steepening with increasing distance from the 
disk. Moreover, one may also argue for a smooth radial 
dependence for the parameter 6. Some theoretical argu¬ 
ments are found in [23, 34]. Observationally, a spectral 
steepening toward high latitudes or longitudes has been 
found in y-rays by the Fermi-hAT observatory [54]. As 
argued in [55] and very recently shown in [56] , the Fermi- 
LAT data are suggestive of a spatial change of CR dif¬ 
fusion. These works are nicely complementary to the 
present paper. 

Leptons — In contrast to hadrons, light CR leptons are 
subjected to synchrotron radiation and inverse Comp¬ 
ton losses with a characteristic time-scale t{E) oc E~^ . 
These processes limit the range that they can travel to 
distance scales A^o ^ \/Tki/o, in the two propagation 
regions. Above a few GeV one has Xo < C in the outer 
halo so that, in contrast to hadrons, CR leptons de¬ 
tected at Earth have spent a larger fraction of time in 
the inner halo. At higher energies the distance scale de¬ 
creases further. The resulting effect is that the spec¬ 
tra above a few tens of GeV are essentially determined 
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FIG. 4. Energy spectra of positrons and electrons multiplied 
by , and positron fraction e"'"/(e“ +e^). The model calcu¬ 
lations are shown in comparison with the data [5, 6, 59-62]. 

by the propagation in the inner halo. Radiative losses 
also steepen the spectrum with a stronger effect than 
diffusion. The lepton fluxes at E >10 GeV are of the 
type (j)e± ~ {t and their spectral index scales as 
7 e ss —Vf. — i(i5 -|- 1) where, for the two-halo scenario, we 
have K = ki and 5 = 5i. In summary, the leptonic spec¬ 
tra under the two-halo model are expected to be harder, 
due to propagation, by a factor i((5 — 5i) Ri0.15-0.2. In 
addition, secondary have a further hardening factor 
coming from the spectra of their parent nuclei p-He. The 
results of our numerical calculations are shown in Fig. 4. 
The force-field approximation is inadequate for CR lep¬ 
tons [57], hence we focus owE ^ GeV. It can be seen that 
our model gives an enhanced flux of secondary positrons 
at high energy. For primary electrons we expect a simi¬ 
lar effect, but their flux depends on the adopted injection 
index which is not known a priori. Using the quantity 
A$g± = (/)g- — (j}f,+ to minimize the secondary and extra 
components [58], we derived i/g =2.6 for the standard 
model and iZg =2.77 for the two-halo model. Hence the 
propagated electron fluxes for the two models are similar. 
It is interesting to note that our model, with the adopted 
parameters, matches the conditions for having a positron 
fraction that increases with energy. 

As for the case of antiprotons, a high-energy enhance¬ 
ment of the positron flux is a characteristic feature of 


the two-halo scenario of GR diffusion. The high-energy 
positron flux is significantly larger than that arising from 
standard models. For dark matter models with TeV scale 
masses, the increased level of background, together with 
a much slower propagation of the signal, will changes the 
intensity and the spectral shape of the expected yield 
from exotic sources. At this point, one may be tempted 
to explain the data entirely with a purely diffusive 
picture. However, we found that even the extreme choice 
Si =0 is unsatisfactory for GR leptons. Besides, this 
would require exceedingly steep injection spectra and a 
worse agreement with the B/C ratio data. Our choice of 
Si = 0.15 represents at present the best agreement be¬ 
tween B/C data and anisotropy requirements. Nonethe¬ 
less, the current B/C data are not very constraining. 
Light-nuclei data from AMS will allow for much more 
stringent predictions. 

Conclusions — Several observed properties of the CR 
spectrum can be understood in terms of spatial changes 
of CR diffusion in the Galaxy, in contrast to stan¬ 
dard propagation models that assume a unique diffusion 
regime everywhere. In our scenario, a departure from the 
standard power-law form is predicted for all elemental 
fluxes. We found an excellent agreement with the new 
AMS proton data at 1 GeV - 2TeV of energy. Our 
model can reconcile the steeply falling low-energy shape 
of the B/C ratio with the high-energy requirements of 
the anisotropy data. In this we have devised perhaps one 
of the few viable solutions for the anisotropy problem. 
Our model shares some consequences with the model pro¬ 
posed in [19]. The y-ray diffuse emission may be useful 
to discriminate the two scenarios. Our predictions have 
also implications for the Galactic synchrotron emission, 
which traces the electron spectrum at different latitudes 
[63]. Interesting results are found for the antiparticle 
spectra, especially in view of recent excitement about 
possible evidence for dark matter. Our model predicts 
remarkably hard spectra for secondary antiproton and 
positron at high energy. In particular, the p/p ratio is 
found to flatten above ^ 10 GeV, while the positron frac¬ 
tion is found to increase with energy. These results have 
a direct impact on upper limits on new physics signals 
using the p/p ratio, or on dark matter models for the 
positron fraction. The antimatter spectra are strictly re¬ 
lated to secondary nuclei. The main limitation is repre¬ 
sented by the large uncertainties of the pre-AMS nuclear 
data. A more robust assessment will be done upon the 
forthcoming release of the AMS data on He Li, C, B/C 
and p/p. In particular, the Li flux will be a powerful tool 
at the level of precision expected by AMS. If new data 
confirms our scenario, then the current standard models 
of GR diffusion will have to be revised to account for it. 
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